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Abstract. We present the results of the planar diffusion of a Dirac particle by step 
and barrier potentials, when the incoming wave impinges at an arbitrary angle with the 
potential. Except for right-angle incidence this process is characterized by the appearance 
of spin flip terms. For the step potential, spin flip occurs for both transmitted and 
reflected waves. However, we find no spin flip in the transmitted barrier result. This is 
surprising because the barrier result may be derived directly from a two-step calculation. 
We demonstrate that the spin flip cancellation indeed occurs for each "particle" (wave 
packet) contribution. 

Subject Classification: 03.65.Pm (pacs). 

I. INTRODUCTION 

In a previous paper [T] , we investigated the one-dimensional phenomena of diffusion of a Dirac particle 
from step and barrier potentials. One of the first observations made in that work was the simplifying 
fact that spin flip did not occur for either potential. This result was not limited to the non relativistic 
limit where it might have been expected. It is however, as we shall show below an exceptional result. 

We consider in this paper potentials and particle rays situated in a plane, the y-z plane. The 
potentials are functions of only one variable, i.e. V = V{z), while the incoming particle direction lies 
in the y-z plane with 9 the impact angle with the potential. The outgoing wave momenta, be they 
reflected, transmitted or those in the barrier region, must of course also lie in the y-z plane. 
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In general, for arbitrary 0, spin flip eontributions oceur for any given incoming polarization. Only 
when = 0, i.e. for right-angle impingement does spin flip completely disappear, and this applies both 
to the step and barrier potentials. This limit case reproduces exactly our previous one-dimensional 
results [T]. 

Our primary objective in this paper is to present the more general planar results, i.e. those for 
arbitrary 0. In doing so, we observe that spin flip is rigourously absent for the transmitted wave in the 
case of the harrier potential. This is not a resonance phenomena, for even in the so called "particle" 
limit, in which infinite transmitted (and reflected) waves occur, no spin flip is found. This limit is 
that pertaining to an incoming wave packet which is small in z compared to the barrier width (L). 
We shall calculate this particle limit by means of the two step method [2l[3] and explicitly verify the 
absence of spin flip in transmission. The sum of the infinite individual waves reproduces the plane 
wave results (maximum interference) including resonance phenomena. 

Before passing in the next section to a detailed discussion of our planar diffusion, we wish to recall 
here that for the one-dimensional Dirac equation three forms of interaction occur with a potential, 
be it a step or barrier of height Va, depending upon the energy of the incoming particle E. For 
E > Vo + m diffusion (2?), ioT Vo — m < E < Vq + m tunneling barrier (T), and ioY E < Vo — m 
Klein energy zone (/C). In the diffusion [1] and Klein [4j[5] cases oscillatory solutions exist everywhere. 
Whereas the tunneling case [SlIT] is characterized by real exponential solutions in the potential region. 
The interpretation of the Klein case conventionally involves pair production as an interpretation of 
the Klein paradox [8]. Only the diffusion case interests us in this paper although most (but not all) 
formulas can be analytically continued into the other energy zones. 

II. GENERAL CONSIDERATIONS 

The time independent Dirac planar {y-z) equation for potential "scattering" is easily reduced to a 
one-dimensional problem when the potential is only a function of one (say z) of the planar variables. 
Separation of variables results in 

^'(r, t) = V(z) cxp[i(p22; - Et)] , 

where is the momentum component along the y-direction and remains so in all regions of the plane. 
In free space, V = Q, V'(^) satisfies 

— i {z) -f Qfaft V'('Z) + /3 m'ip{z) = E Tpi^) ■ (1) 

With ps the z-componcnt of momentum, we have as one of the explicit solutions (polarized in the 
z-direction) 

ipiz) ^ [1, 0, ps/iE + m), ip^/{E + m)Y exp[zp3z] , 

up to an overall normalization factor. We note that p^/ps = iaxiO and E = \/p\ + p\ + rn^ ■ Through- 
out this paper we consider the incoming particle (travelling from negative z) polarized as above, but 
our results will be indifferent to the specific choice of polarization and indeed will be expressed in 
terms of spin conserving and spin fiip. 

For the region in which V{z) — Vo, we must make the following modifications (translation) 

E ^ E -Va , P2 P2 , Pi^ , 

with {E — Va) = P2 + 93 + ''^^^ ■ The existence of a non zero value for p^ modifies our spinors when 
compared to our previous one-dimensional calculations, and hence our diffusion results It also 
complicates the kinematic conditions for being in one of the zones 2?, T and /C. The simplest way to 
see this is to define a new mass , 

= ^pI + = \/ E^ sin^ + cos^ 9 . (2) 

Then, the one-dimensional energy zones are simply generalized to 

(2?) : E>Va + , 

(T) : Va~m^ < E <Va + m^ , 
{iq : Va~m^> E . 
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In Figs. (1) and (2), we plot the separation of the various zones by fixing respectively the incident 
angle and the potential Vq. The Klein zone is absent for VJ, < 2m, see Fig. (2). The plots in Fig. (2) 
also show that for a given energy E, we may transit for high E from 2? to T or for low E from /C to 
T by varying the incidence angle 9. For Vo — m < E < Vo + only the T zone exists independent 
of 0. We note that by varying 9, we can never pass through all three zones. We also note that the 
one-dimensional kinematics are simply given by the energy points on the axis sin^ = in Fig. (2). 

The Dirac equation, being of first order in the spatial derivatives, implies that continuity equations 
are applied only to the field. However, this provides four equations, one for each component of the 
spinors involved. The "small" components correspond in the non relativistic (NR) limit to the continu- 
ity of the spatial derivative of the NR field. The remaining doubling of continuity equations (compared 
to the Schrodinger equation) is exactly what is needed to determine both the spin conserving and 
spin flip contributions. 

In the next section, we apply the continuity condition to the step at z = for waves coming from 
either direction (needed in the two-step procedure) and at 
the step calculations of planar diffusion: 
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In these figures, we list the contributing reflection (i?) and transmission (T) amplitudes. The ± 
subscript distinguish an up step from a down step while the tilde represents reflection and transmission 
for an incoming wave from the right. The prime terms correspond to spin flip amplitudes. In the next 
section, we also give the results of the continuity equations for the barrier. 



III. PLANE WAVE RESULTS 

Consider the step potential defined by 

Viz) = { for z < , Vo for z > } , 

with Vq > and with an incoming plane wave from the left (see the step A) with a definite polarization 
(spin along the z-direction) . The spinor continuity equations are 




exp[«p3z;J 



E^ 



i?+(0) 



E 



T+(0) 




cxp[- ipsZj 



CXp[ IQsZ \ 



(3) 



These matrix equations can be rewritten as 



{E + m)[l + R+{0)] 
(E + m) R'+{Q) 

Ps[l - R+m - iP2R'+m 
ip,[l + R+iO)] + p,R'^{0) 



E + m 

E-V„ + ', 



{E-V„ + m)T+{0) 
{E-Vo+m)T[{0) 
q,T+{0) - iq,T[{0) 
iiq,T+{0)^q,T!^{0) 
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Solution of which yields 



i?+(0) 
T+(0) 
i?',(0) 

t;(o) 



{pI + to^ + mE) K 



P2P3 Vo 



{E + m){pl 



VoE) 



VoE) ' 



(4) 



Here, we observe that because the step is situated at z = 0, the i?+(0) and T+(0) amplitudes are 
real while the -R+(0) and T'^(O) are imaginary. All four amplitudes exist. If we change the incoming 
polarization (from spin up to spin down), we find exactly the same solutions, although from a different 
form of matrix equation. This is a general property and we shall henceforth not repeat this observation. 

As a simple check of our results, we note that when K only T+(0) = 1 survives as must be. We 
also observe that as — >■ {9 ^ 0) the spin flip terms vanish and this confirms our one-dimensional 
results. The other two step cases (B and C) can be calculated in the same way and yield the results 



R-{L) 
T-{L) 
R'_{L) 



[pl + TO^ + m{E - K)] K 
[E -Vo+m){pl+p,q,-VoE) 
ql{E + m)+ p,q,{E -~ Vo + m) 



m){pl+P3qs 



VoE) 



[E^Vo 
R'_{L) cxp[- 



exppiqai] , 
exp[i(g3 -P3)L\ , 
cxp[2i93L] , 



- m){pl +psq3 - VoE) 
-i(q.,+Ps)L] = R'_{0) cxp[i(g, 



T'_{L) 

While for right impingement at z = wc obtain, 

^_(0) = i?_(0) , ^^(O) = r_(0) , ^'_(0) = -i?'_(0) , f'_{Q) = R'_{Q) 



(5) 



(6) 



Passing now to the barrier, we define the reflection and transmission amplitudes respectively in 
regions I and III by i?, i?', T and T', while for region II, that of the potential Vq, we use A, A', B 
and B' . These are not related in a simple manner to our previous step results. The connection will 
however be derived in the next section. 

The expression for ip{z) within each region are given below. Region I (z < 0): 



f E + m 
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Pi 

iP'i 



cxp[ip3Z 




cxp[ — ip^z \ 



(7) 



Region II (0 < z < L): 



Region III (z > L) 
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(8) 
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In the above expressions factors such as !/(£' + m) and 1/{E — Vo + m) have been absorbed into 
the amphtudes for simphfication. After ehmination of the intermediate A, A' , B and B' the coupled 
continuity equations yield the matrix equation 



{E + m){l + R) ■ 

{E + m) R' 
P,{1 - R) - ip,R' 
ip^il + R) + p,R' 



M 



■ {E + m)T ' 
{E + m) T' 
PiT - ip2T' 
ip2T - psT' 



exp[ip3L] 



where 



M = — 

93 



-i{E 



93 cos(g3L) 
92 sin(g3L) 
-Va-m) sin(g3L) 

i{E 



ga sin(g3L) 
93 cosiq^L) 


-Vo - m) sm{q3L) 



-i{E — Vo + m) sin(93L) 

i{E 
93 cos(93L) 
92 siu(93L) 





-Vo + m) sin(93L) 

92 sin(93L) 

93 cos(93L) 



The solution of these equations are 




P2 



Vo 



sin(93L) 



E + m J ■ " PsqsJ^ 
smiqsL) 



T 



E + m ■ " P393J^ 
cxp(— zp3L) 



T' = 



(10) 



with 



. P3 - EVo 



sin(93L) 



These are the generalized barrier results for a plane wave. They contain the momentum p2 indicating 
a dependence upon incident angle. For p2 = 0, we reproduce the one-dimensional Dirac barrier results 
published in our previous paper [1]. Spin flip is indeed absent in this limit. However, the surprising 
feature of the general barrier results is that T' = for all incident angles. This was not expected and 
seems to contrast with the fact that the equivalent terms T!j_ are not identically null for the steps. 

In the next section, we perform the two step calculation to redetermine the above expressions, in 
particular that for T' . This will confirm the above results and demonstrate that T' = is not a reso- 
nance phenomena. On the other hand typical resonance behavior is present in the above expressions. 
Whenever sin(93L) = 0, i.e. for q^L = wk (with n a positive integer) both reflected amplitudes, R 
and i?', vanish and the transmitted probability \T\ = 1. Plane waves are theoretical abstractions, 
the above results are in truth good approximations only for harrier widths much smaller than the in- 
coming wave packets widths. In other words high values of n corresponding to large L will not exhibit 
resonance behavior [T]. 



IV. TWO STEP CALCULATION 

In this section, we recalculate the barrier results using the step results. Specifically, this approach for 
the barrier is called the two step method and does not directly involve the continuity equations of 
the previous section. The method uses three step results. In addition to that for the step ai z = L 
impinged upon from the left, it also uses the step results at z = twice. Once for the initial incoming 
wave (impingement from the left) and then those for a wave reflected from the end of the barrier with 
consequent impingement at z = from the right. 

This method of calculation can be used to reproduce the standard barrier results by simply adding 
the infinite contributions to transmission and reflection yielding the so called (plane) wave limit. By 
treating each contribution as incoherent with the others we obtain the particle limit. Probabilities are 
conserved in both limits although the total transmission/reflection probabilities are quite different. 
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The wave limit for example is characterized by resonance phenomena, the particle limit is not. In 
addition to calculating the individual particle limit probabilities for transmission, it is one of our 
objectives to control if vanishing spin flip is a resonance phenomena or not. 

The incoming wave is polarized. At z = it encounters the first of the two steps and two con- 
tributions proceed (are transmitted) to the right. These are indicated by the spin conserving Tj^{Q) 
and the spin flip term T'^(O). At the second (downward) step at z = L each of these contributions 
produces two transmitted amplitudes. These four transmitted terms combine into two sums: the spin 
conserving amplitude, 

A, = T+(o)r_(i) + r;(o)r:(L) , (ii) 

and a spin flip amplitude, 

a; = T+(o)r(i) + r;(o)r_(L) . (12) 

Now, we can use the step results previously given to observe that 

A^= '^P''^' exp[»(g3-P3)X] and ^',=0. (13) 

PI + PsQa - EVo 

So the first particle spin flip contribution is null. This fact alone docs not guarantee that subsequent 
spin flip contributions are null. For example, the second transmitted contributions contain in addition 
to the T factors also two R factors corresponding to an additional back and forth passage over the 
barrier, 

A, = i?_ (L)i?_ (0) + {L)R'_ (0) = exp[2 ^q,L] , (14) 

A'^ = R- {L)R'_ (0) + R'_ {L)R^ (0) = . (15) 
The overall second contribution to the spin flip is thus 

^T^n A^An = . 

All higher spin flip contributions take the form of the the above second contribution multiplied by 
powers of the spin- conserving double reflection factor Thus, the vanishing of the second spin flip 
term does indeed imply the vanishing of all the spin flip contributions. 

We now calculate the individual spin conserving contributions. The first contribution A^. has 
already been given (1111) . The second contribution is 

All higher (later emerging) contributions are now obvious. The n-th term reads, 

A-tA^ 

In the wave limit these contributions must be added coherently to give a single outgoing transmission 
amplitude. This may conveniently be written as, 

(16) 



1-^R 



After inserting the specific step expressions of the previous section, we reproduce after a little algebra 
the plane wave result for the barrier transmission, 

y ^ ^ cxp(-zp3£) 
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V. CONCLUSIONS 

In the literature one normally encounters one dimensional potential analysis. However, when spin 
and relativity are relevant one-dimensional analysis may be too limited. For example the absence of 
spin flip terms for one-dimensional step and barrier diffusion is no longer valid for planar diffusion 
where an angle of incidence exists. The potentials are still considered functions of a single spatial 
variable (z in this paper) but the incoming particle have two momentum components. This small 
modification produces significant differences, specifically the appearance, in general, of spin flip terms 
and consequent modifications in the non-flip amplitudes. We have demonstrated these facts in this 
paper. These effects are a direct consequence of the angular dependence of the Dirac spinors. 

We have found one notable exception to the above, for which we have no simple explanation, 
although we believe one must surely exist. This exception is that in the case of the barrier potential 
there is no spin-flip transmission amplitude. There are always spin-flip terms for the step, be the step 
rising or dropping , except of course for the one dimensional limit when the potential is met head on. 
This makes the barrier result all the more unexpected, since we know that the barrier result may be 
derived from a double-step analysis which uses only the step results. 

We have listed in the previous sections all relevant planar diffusion results for the step and barrier. 
We have also derived the transmission amplitude for the barrier in the " particle" limit via the two-step 
method. This demonstrates that the absence of spin-flip for each individual outgoing wave packet, 
independent of the degree of coherence involved, i.e. it is not a resonance type phenomena. 

As a side product, we have described the kinematics of dirac planar scattering and observed that 
by varying the incidence angle, for a given incoming energy we may transit through two kinematic 
zones, e.g. from diffusion to or from tunnelling, but never through all three kinematic zones P, T, /C. 
We have in the past referred to these zones as energy zones, but this is correct only for one-dimensional 
studies. They should instead be referred to as kinematic zones which depend upon both energy and 
angle of incidence. Energy alone, normally does not determine these zones. These kinematic zones are 
distinguished by quite distinct physics. In this paper, we have limited our attention to diffusion. 
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Fig. 1. The separate kinematic zones for three choices of incoming angle as a function of potential height. 
Case (c) corresponds to the one-dimensional plot. 
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Fig. 2. The angular dependence separated into kinematic zones for three choices of potential. Case (c) is the 
limit case (highest potential) for which there is no Klein zone. 



